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Using an advanced version of the renormalized mean-field theory (RMFT) for the t-J model, we 
examine spin-singlet superconducting (SC) state of d^2_y2-symmetry. Overall doping dependence 
of the SC gap magnitude is in good agreement with experimental results for Bi2Sr2CaCu2 08+i 
(BSCCO) and La2-a:Sra;Cu04 (LSCO) compounds at the optimal doping and in the overdoped 
regime. We also calculate the dispersion relation for the Bogoliubov quasiparticles and compare 
our findings both with the angle resolved photoemission data for the cuprates, as well as with the 
variational Monte Carlo and other mean-field studies. Within the method proposed by Fukushima 
[cf. Phys. Rev. B 78, 115105 (2008)], we analyze different forms of the t-J Hamiltonian, i.e. 
modifications caused by the form of exchange interaction, and by the presence of three-site terms. 
It is shown that although the former has a small influence, the latter suppresses strongly the su- 
perconductivity. We also analyze the temperature dependence of the gap magnitude and compare 
the results with those of the recently introduced finite-temperature renormalized mean-fleld theory 
(TRMFT) of Wang et al. [cf. Phys. Rev. B 82, 125105 (2010)]. 
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I. INTRODUCTION 

One of the most characteristic features of high temper- 
ature superconductivity (SC) is that upon hole doping, 
with the hole concentration x > 0.05, a generic antiferro- 
magnetic Mott insulating state of e.g. La2Cu04 -'^ trans- 
forms into SC state. The latter, in turn, after reaching a 
maximal transition temperature at a; 0.15, disappears 
at the upper critical concentration Xc ^ 0.25 — 0.35, de- 
pending on the systemj^id The last property is particu- 
larly surprising, since in the overdoped regime x > 0.15 
the system evolves gradually from a non-Fermi liquid 
into a quantum liquid that can be regarded as an un- 
conventional Fermi liquid.^ The appearance of Xc may 
speak in favor of real-space type of pairing, as the in- 
creased hole doping reduces the pairing correlations in 
real space. To describe the above features the t-J model 
is often invokedr'"— and the kinetic exchange interaction 
is claimed to induce both antiferromagnetism and the 
superconductivity. ^ 

The basic question is whether within the renormal- 
ized mean-field theory (RMFT)i^iii we can reproduce 
at least some of the above properties in a semiquanti- 
tative manner,—"— as advocated strongly i^^i^^ On the 
other hand, variational Monte Carlo (VMC) methods, 
within which one treats exactly the double occupancy 
exclusion, is knowr>i^"— to provide a semiquantitative 
description of the SC correlated state. Hence it is of- 
ten regarded as being superior to any mean-field (MF) 
treatment. However, a proper MF approach would have 
important advantages over VMC. First, its results are not 
limited to small clusters. Second, it can offer an analytic 
insight into the physical contents of the model. Third, 
it allows for a detailed comparison with experiment or, 
strictly speaking, its critical assessment, as detailed be- 
low. 



Numerous attempts to improve RMFT have been 
madef22r"22, in order to take into account also intersite 
correlations. Here we show that within the MF renor- 
malization scheme proposed recently by Fukushima;^ 
supplemented with the maximum-entropy based^^ self- 
consistent variational approach^"— we can produce, 
among others, the results that are competitive to those of 
VMC. Specifically, the upper critical concentration 
well as the principal features of the excitation spectrum, 
are shown to agree quite well with experiment in the over- 
doped regime. Finally, motivated by a recent paper 
we comment on the extension of the present treatment 
to finite temperatures. We study a behavior of the renor- 
malized gap magnitude as a function of the temperature, 
which is shown to be in quite good agreement with the 
classic result of the BCS theory. 



II. MODEL AND METHOD 

We start from t-J model^"— is expressed by the follow- 
ing Hamiltonian 

iitj ^P^Ht + Y, (s, • Sj - jhi^,) + C2H3) p. (1) 

The first term, Ht = X^ijo- ^ijcLcjo- is the kinetic en- 
ergy part, the second expresses the kinetic exchange, and 
the third the three-site terms. means the summa- 

tion pair of sites (bonds). The Gutzwiller projector 
P — J^j(l — fii-fhi^) eliminates double occupancies in real 
space. Also, explicitly 
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with the projected fermion operators defined as: bia = 
(1 - niij)cia, i^ia = (1 - n,g)n,cr, = J2a'^i<^' ^iid 
5'°' = bl^big. Here the standard fermion creation (anni- 
hilation) operators are cj^. (cia) and hi^ = c\^Cia. Hamil- 
tonian ([T]) in its complete form, i.e. with ci = C2 = 1 is 
derived by applying canonical transformation to the Hub- 
bard Hamiltonian, Htu, in the strong couphng (|t| << U) 
limit,— 1^ Htj = exp(— iS')_fft(7 exp(iS'). Then, kinetic ex- 
change integral — 4:\tij\'^/U. On the other hand, the 
ci = C2 = case has been used as a simplified effective 
model describing low-energy sector of a multi-band Hub- 
bard model.— 12 Both the latter form, or that with ci = 1, 
C2 = 0, are often used in studies of high- Tc supercon- 
ductivity. Our aim here is, among others, to compare 
the complete version of the t-J model with the simpler 
forms. 

It has been arguedfi^ that a correct variational state 
describing the ground state of the t-J Hamiltonian is 
of the form |^) = P|4'o), where \'^o) is an uncorre- 
lated state, here taken as the Bardeen-Cooper-Schrieffer 
(BCS)-type state. When the Hubbard model is the start- 
ing point of analysis, the re-transformed back state with 
perturbatively reintroduced double occupancies is used 
instead, i.e |*) |§) = exp J^^i^ Also, the de- 

tailed form of l^'o) may be postulated without speci- 
fying the underlying microscopic Hamiltonianji^ How- 
ever, within RMFT, l^fo) may also be selected as an 
eigenstate of the effective single-particle Hamiltonian, 
obtained from the mean-field treatment of ([1]) fiSiiii22ii2i 
This point of view is also taken up here. 

Construction of RMFT requires a prescription (here 
called the renormalization scheme, RS) for calculating 
expectation value of an arbitrary operator O in the cor- 
related state 14"). Explicitly, 



~ _ (*ia|vi/) 



^o\PcOPc\^o) 
(*o|P^|*o> 



{PcOPc 



(3) 

In the above, the projector P has been replaced by 
more general correlator Pc, which differs from P by the 
presence of fugacity factors, ensuring that the projected 
and preprojcctcd average particle numbers are equal, i.e. 

We adopt the RS of Ref. 22 and supplement it with 
the variational formalism proposed recently by us.— 
Combination of these two factors yields the results that 
differ considerably as compared to those of the standard 
RMFTiiSiii To carry out the self-consistent variational 
approach, we start from the following effective MF grand 
Hamiltonian 



K\ = {firja{4a(^j^ - Xrja) + H.C. 

- {"^'i ^ ) + ^■^) ^ XI ^^"^ 

(ij) ia 

- 5](A:r,(n». +VK(x^Ja,A,J,n,,), (4) 



with Kx ^ Hx- ijlN, W = (-fftj)c, = {hia), Xij = 
(c1o-Cj<t), Ay EE (c,tCji-Cj^Cjt)/2, and A,j = (Ay). The 
averages (. . .) are defined by ([3]) for zero temperature, 
or in general, with the help of density operator px, (see 
below), fjijcr, ^ij, and A^, (not to be confused with the 
fugacity factors of Ref. |23 ) are the Lagrange multipliers, 
ensuring the self-consistency of this variational approach, 
as this is not guaranteed within the standard Gutzwiller 
approximation.— The form Q , apparently different from 
the usual formulation of RMFT, is fully equivalent to it 
(a similar approach is discussed in Ref. [23l ). 

We solve this model on a square lattice and in the spa- 
tially homogeneous situation, with no coexisting mag- 
netic order. The model parameters are ti = t,t3 = t' 
and t<^ = i", where s = 1,3,5 corresponds to sites located 
at the distances d{i,j) of 1, ^/2, and 2 lattice constants, 
respectively. Consequently, the following Lagrange mul- 
tipliers and the corresponding mean fields are assumed 
as nonzero: Xij = 6 and 77^- = 77^; = A^(^y-), with 
Ax = A = ~Ay and 7^ = j^^y^ = ±7, both for 

= 1; "■io- = "-iff = n/2 and A"^ = A. Also, in 
^ we retain all the terms of the orders of /U and 
tt' /U, and neglect a smaller term of the order [t'Y /U . 
Then, {Htj)c is obtained using the formalism of Ref. [22I . 
Explicitly, 

W = {Ht)c + {Hj)c + {H^)c = Wt + Wj + W^, (5) 

where Wj = Wj — AJcin^/2, and A is the number of 
lattice sites. In effect, we have that 



A 



2 - n 



s=3,5 



(2 - n)- 



(2 - nf 



(6) 



Wj _ , J 3 (A^ + e,) + ci(l - nf (A^ - \ 

. (7) 

The formula for 14^3, containing three-site terms, is too 
lengthy to be reproduced here. 

The next step is the diagonalization of Hx (SI) via 
Bogoliubov-Valatin transformation, which leads to 

Kx = Y. ^k(7ko7ko + T^iTki) + - i?k) + C, (8) 



with C = iy-HA(8X;s6'7s + 4A7 + An). Also, the quasi- 
particle energy Sk, the renormalized gap I?k, and the 
renormalized band energy are respectively given by 



E^. = \Sl + Dl i?k = -:^r_(k), (9) 



a = - (77ir+(k) + ?73e(k) + r;5r5(k) -i- , (lo) 
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r±(k) = 2(cos(fca;) ±cos{ky)), 0(k) = A cos{kx) cos{ky) , 
r^Ck) = 2{cos{2kx)±cos{2ky)), and /i = ^ + A. Then, the 
corresponding grand-canonical density operator reads 

px = Z-'cM-PKx), (11) 

with Zx ~ Ti-[exp{—f3Kx)], (3 — l/ksT, and the general- 
ized Landau functional, J- = —l3~^ InZ^ reads 

= C + 5]((a-i?k)-|ln(l-l-e-^^'')).(12) 

Because of extra constraints in (U) , the necessary condi- 
tions for J-^ to have a minimum are 

d^T = 0, d,T = 0, (13) 

for the mean-fields w e {^s,A, n}, n = 1 — x, and 
the Lagrange multipliers z g {?7s,7. A}. In contrast to 
our previous formulation)^^ the equations dwJ^ — are 
solved analytically, but the solution of those remaining 
{dzT = 0) must be determined numerically. Optimal val- 
ues of the mean-fields and the Lagrange multipliers, given 
by a solution of characterized with the lowest value 
of T and denoted respectively {^[^\ A'^^\ n) = 

MT,V,fi), iv?\vi°\v^°\l^°\>^^"^) ^ MT,V,fi), may 
be used next to construct standard grand thermodynamic 
potential ft and the free energy F = + . Explicitly, 

n{T, V, n) = T{T, V, 11- io(T, V, ii)MT, y))- (14) 

The present formalism, based on the maximum-entropy 
principle^** is formally valid for arbitrary T > 0. Conse- 
quently, we have replaced pure state jV'o) by mixed state 
represented by p\. However, the solutions obtained for 
non-zero, but sufficiently low T, are for all practical pur- 
poses identical to those for T = 0. This situation is stud- 
ied in the next Section. The separate question is that of 
the validity of RMFT approach at T > from the point 
of view of physical consistency. This is analyzed in more 
detail in Section IV. 



III. RESULTS FOR T 

An analytical expression for W = {Htj)c allows us to 
make qualitative predictions before the numerical anal- 
ysis is carried out. First, from ([5]) we expect a strong 
tendency to the superconductivity suppression for higher 
doping, as SC order leads to the band energy decrease 
~ A'^. On the other hand, in that regime the renor- 
malized band energy becomes predominant over the ex- 
change part. In effect, the normal state is favored over 
SC for X > Xc with Xc smaller than obtained within pre- 
vious MF treatments.— "—ii^ Second, from (O we infer, 
that the influence of the ViVj/A term on Wj is small ex- 
cept for the largest doping; this is due to the presence of 
the (1 — n)^ prefactor (the other term ~ merely shifts 



the chemical potential). On the other hand, W3 is mul- 
tiplied only by (1 — n) = x prefactor. Consequently, for 
higher x this term becomes rather important, due to the 
number of distinct three-site terms present for a given 
initial site and spin direction (eight for d{i,j) = and 
four for d{i,j) = 2). Also, this part of W is expected 
to suppress SC order, as the term ~ in W3 contains 
a factor, which is positive for reasonable values of other 
mean fields and model parameters. 

We solve numerically the part dzJ- = oi Eqs. 
for the mean fields using periodic boundary conditions on 
the lattice of A = 512^ sites, to minimize finite size effects 
and for low temperature fc^T = 2 • 10~^J. The solution 
amounts to solving simultaneously the system of five non- 
linear equations using GNU Scientific Library (GSL). In 
the most cases we take the parameters \t\/J = 3 (corre- 
sponding to U/\t\ = 12 for the Hubbard model), t'/t — 
or (—0.25), and t" = 0. Additionally, we take also 
\t\ — 0.3 eV or \t\ — 0.4 eV, which correspond roughly to 
the lower and the upper limits of the realistic values of 
this parameter, depending on the compound. Values of 
\t\ close to 0.4 eV have been determined from the band- 
structure calculations;^!^ whereas |t| = 0.3 eV is used in 
Refs. mill To highlight the infiuence of various forms 
of ([T]), the results for different values of ci and C2 are 
analyzed. The numbers 1, 2, 3 (4, 5, 6) in Figs. [T]|4]and 
Table I correspond to the three situations: ci — C2 ~ 
(i.e. with the Sj • Sj part of the kinetic exchange only, 
cf. H])), ci = 1 and C2 — (with full form of the kinetic 
exchange), and ci = C2 = 1 (complete form of the t- 
J Hamiltonian with the three-site terms included), each 
case taken for t'/t = {t'/t = —0.25), respectively. 

In Table I we provide the equilibrium values of the 
mean-fields and of the Lagrange multipliers for cases 1-6 
and for x = 0.175, a representative hole concentration in 
the overdoped regime. 

Table I. Optimal values of mean-field parameters for 







T w 


and X 


= 0.175 








1 


2 


3 


4 


5 


6 




0.1970 


0.1969 


0.1990 


0.1924 


0.1922 


0.1944 




0.0468 


0.0465 


0.0505 


0.0241 


0.0239 


0.0225 




-0.0080 


-0.0076 


-0.0144 


0.0337 


0.0340 


0.0383 


A, 


0.0687 


0.0708 


0.0202 


0.0903 


0.0919 


0.0534 


m 


1.0080 


1.0030 


1.2355 


1.0031 


0.9982 


1.1845 


m 


0.0000 


0.0000 


0.0803 


-0.2223 


-0.2223 


-0.2118 


V5 


0.0000 


0.0000 


0.0408 


0.0000 


0.0000 


0.0404 




0.1584 


0.1665 


0.0320 


0.2126 


0.2205 


0.0834 


ft' 


-0.4069 


-0.4080 


-0.2935 


-0.8633 


-0.8614 


-0.9406 



In Fig. [T] we plot the dispersion relation for the Bogoli- 
ubov quasiparticles, calculated for the parameters dis- 
played in Table I. The infiuence of on i?k is of com- 
parable magnitude to that of having nonzero t' . 

Next, we discuss the doping-dependence of the renor- 
malized SC order parameter {^ij)c = (cf. Eqn. (18) 
of Ref. [2^ . The numerical results confirm the above 
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FIG. 1: (Color online) Dispersion relations along the main 
symmetry lines in the two-dimensional Brillouin zone for 
n = 0.825. The various curves are explained in main text. 
Note, the curves 1 and 2, as well as 4 and 5 are practically 
indistinguishable (the influence of the term '--^ ci in ([1]) is 
neglible). 

made qualitative predictions. In Fig. [2] we plot Ac for 
the cases 1-6 specified above, as well as for t'/t = —0.27 
(value being reasonable for BSCCO compoundSf22) and 
J/\t\ — 0.3 (curve 7). Note, that the upper criti- 
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FIG. 2: (Color online) Doping dependence of the renormalized 
superconducting order parameter (A)c = Ac. The curves 1- 
6 correspond to those in Fig. 1. The curve 7 is for ci = 1, 
C2 = 0, t'/t = -0.27, t" = 0, and for J/\t\ = 0.3. For 
definitions of cases 1-7 see the main text. Cases 3 and 6 
correspond to the complete t-J Hamiltonian ([1} with ci = C2. 

cal concentration Xc for the cases 4 and 5 is close to 
the VMC result,^'' obtained within the Hubbard model 
and using \ip) for the corresponding model parameters. 
Also, nonzero value of t' enhances superconductivity, 
in agreement with previous VMC results^^ and other 
calculations 1^ Let us emphasize again, that the presence 
of the iJa term acts in the opposite direction. The verti- 
cal line roughly marks the boundary between under- and 
over-doped regimes. Importantly, for \t\/J = 3, t" — 
and different t'/t values, < t' /t < 0.25, Xc lies in the in- 
terval 0.2 < < 0.35, depending on the form of Htj, as 
illustrated in Fig. [5] As said above, the small difference 



between the curves 4 and 5 (as well as between 1 and 2) 
shows an insignificane of the term ~ ci (cf. also Table I). 
The results for Xc are in a good overall agreement with 
the experimental data for the cupratesj^'S, To provide an 
additional support for our results, we list in Table II the 
values of function of J, for either t'/t — —0.1 

(considered to be relevant to the LSCO compound^), or 
t'/t = -0.27. 

A remark is in place here. The family of curves in Fig. 
[H and of Xc values in Table II has the following meaning. 
Each Xc value singles out either the choice of a model or a 
particular set of parameters. This detailed analysis is to 
illustrate that there is a clear upper critical concentration 
in the proper range, irrespectively of the model details or 
particular set of parameter values. 

Table II. Upper critical concentration Xc vs. J for 
t' = -0.27i or t' = -O.lt. The symbol A (B) labels the 
case C2 — and ci = (ci = 1), respectively, whereas C 
means that ci = C2 = 1 is taken in the computation. 



J/\t\ 


0.2 0.3 0.333 0.375 0.4 


t'/t = -0.1 A 
t'/t = -0.27 A 


0.18 0.26 0.29 0.32 0.33 
0.2 0.31 0.34 0.38 0.4 


t'/t = -0.1 B 
t'/t = -0.27 B 


0.18 0.27 0.3 0.33 0.35 
0.2 0.33 0.36 0.4 0.42 


t'/t = -0.1 C 
t'/t = -0.27 C 


0.15 0.21 0.22 0.24 0.25 
0.15 0.23 0.26 0.28 0.29 



In Fig. [3] we plot x-dependence of the SC gap Dk for 
k = (tt, 0) (cf. Eqn. ^) and compare our results with 
the experimental data.^^i For the selections of t'/t and 
J/\t\ as in Fig. [U no fully satisfactory agreement with 
experiment is achieved in the entire range of x. However, 
the agreement with experiment is quite good for the pa- 
rameters corresponding to the curves 1, 2, 4 and 5 in the 
overdoped regime, both for |t| = 0.3 eV and \t\ — 0.4 eV. 
The best overall fit is achieved for the set of parameters 
represented by curve 7. Note, that in all the cases 1-7 the 
quasiparticle energies obtained here are decisively lower 
than those in the standard RMFT formulationi^^ These 
differences are caused by both the particular selection 
of the renormalization scheme, as well as by the varia- 
tional method we use. As a consequence, we obtain also 
lower values of the Fermi velocity, vp = |Vfe^k||k|=fci7j 
calculated for the nodal ((0,0) — > (7r,7r)) direction. The 
lattice constant has been taken as ao = 

The x-dependence of vp is detailed in Fig. 2] for the 
same set of parameters as in Figs. [5] and [31 for both 
\t\ = 0.3 eV and \t\ = 0.4 eV, and compared with the 
data discussed beforeJ^ The theoretical values are still 
too low. Also, the x-dependence of both -Dk=(Tr,o) 
Vp, obtained within the MF approaches, is stronger than 
observed in experiment. This feature is shared with the 
other mean-field approaches .^^^^ However, the experi- 
mental values for BSCCO 1.5 eV A have also been 
reported," and are quite close to our results. Note that 
the disagreement is largest in the underdoped regime (to 
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directly for T > 0. It may seem, that having a finite- 
temperature formaUsm at our disposal, we may exam- 
ine arbitrary temperature by simply changing value of /3 
in J^(/3) Eq. Unfortunately, this may lead to the 

situation similar to that encountered in the slave-boson 
mean-field theories, which application at finite tempera- 
tures is invalidated by incorrect evaluation of the entropy 
part of the free energyi^ 

Recently, an attempt to extend RMFT to T > have 
been madei^ Within this approach, termed finite tem- 
perature RMFT (TRMFT), the term 



.12 0.16 0.2 0.24 0.08 0.12 0.16 0.2 0.24 
X = 1-n X = 1-n 



AS* = — ( ej In — + Qi In — + di In 

^ ^ eio Qio dio 



(15) 



FIG. 3: (Color online) Doping dependences of the SC gap 
Dk at k = (tt, 0) for cases 1-6 and for t' /t — —0.27, and 
J/\t\ — 0.3 (filled diamonds). Large filled circles - experi- 
mental data.— Note, that in contradistinction to Ref.— no ad 
hoc introduced scaling factor a = 1/2 is necessary to obtain 
a reasonable agreement in the overdoped regime, i.e. to the 
right of the vertical line. 



the left of the vertical line). 




0.1 0.2 0.3 
X = 1-n 



FIG. 4: (Color online) Doping dependence of Fermi velocity in 
the nodal ((0, 0) (n, n)) direction. Experimental data (cf.^"^ 
and References therein) are marked by diamonds (YBCO), 
squares (LSCO) and solid circles (BSCCO). 



IV. 



RESULTS FOR T > 0. 



So far, we have analyzed only low temperature (3 J = 
UbT/J = 500, practically equivalent to the true T = 
situation. Obviously, it would be interesting to extend 
the analysis for higher temperatures, and in particu- 
lar, to determine the critical temperature Tc as a func- 
tion of doping. Standard RMFT approach, as based on 
the Gutzwiller approximation (GA), is devised to study 
ground-state properties and as such, is not applicable 



is added to the single-particle entropy 5*0 = — Trpolnpo 
of the mean- field model (here, po = p\). In Ec^. (|15|) . 
we have = {Ei)c, e.,o = {Ei), = {Qi)c, Qm = {Qi), 
di = (A)c, and d^o = (A) with Ei = {1 - ni-f){l - fii^), 
Qi = fii\{^ - nil) + ~ ^it) and A = ni^nii- Note, 
that AS* < 0. Derivation of (IT5|) . and its possible gen- 
eralizations will be discussed elsewhere.^'* Here we sim- 
ply adapt Eq. P3|) . which may be treated as a reason- 
able Ansatz, similar in spirit to the finite-temperature 
extensions of the Gutzwiller approximation proposed 
earlier ^""i22 

Within the present formalism, AS can be included by 
replacing IF(^5, A,n) (O by W-TAS. However, for the 
t-J model, di — 0, and consequently, for non-magnetic, 
homogeneous solutions studied here, AS = AStj de- 
pends only on the total particle number n, i.e. 



ASt., 
A 



= (2 - n) In (l - I) - (1 - n) ln(l - n). (16) 



Therefore, the presence of AStj results only in different 
values of p. A, and thermodynamic potentials V, and F 
([T4|. Still, fL = fi + X, all the mean fields, and the re- 
maining Lagrange multipliers, as well as the free energy 
difference between the superconducting and the normal 
solutions, remain unchanged. Hence, AS in the form 
(|15p does not lead to any nontrivial modifications of our 
original formulation. 

Nonetheless, it still seems to be interesting to apply the 
above formalism to examine nonzero temperature situa- 
tion. This would reveal limitations of the present form of 
RMFT and should help in formulating more satisfactory 
finite-temperature mean-field treatment of the t-J model. 
First, in Fig. [S] we plot temperature dependence of 
(renormalized) gap magnitude Ac (T) as a function of T 
for selected hole concentrations. For the values of model 
parameters used here, we have J — lOOmeV « 1160K 
(J = 133meV « 1550K) for t = 300meV (t = 400meV), 
respectively. Therefore, the critical temperature is over- 
estimated by a factor 3-5 (e.g. Tc — 340 — 450K at the 
optimal doping and Tc ~ 25 — 35K at a; = 0.3. This 
is a common feature of all mean-field type approaches, 
clearly caused by neglecting the fluctuations. This also 
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FIG. 5: (Color online) Temperature dependence of renormal- 
ized superconducting order parameter (A)c(T) = Ac{T). 
The numbers 1-8 correspond to values of doping equal 0.125 
(optimal doping, o.d.), 0.15, 0.175, 0.2, 0.225, 0.25, 0.275 and 
0.3, all for ci = 1, C2 = 0, t'/t = -0.25. 

shows the insufficiency of the expression AS in the 
present t-J model case. 

Nonetheless, one interesting property of the present ap- 
proach should be noted. Namely, an uncorrelated wave 
function I'I'o), or a related density operator p\ (jlip has 
an essentially the same form as that coming from the 
Bardeen-Cooper-Schrieffer (BCS) theory. Therefore, it 
seems natural to compare our values of reduced renor- 
mahzed gap magnitude Ac(T)/Ac(0) with the standard 
BCS result^ given by A(T)/A(0) = tanh(A(r)/tA(0)), 
where t = T/Tc- The results are shown in Fig. |6l 
Interestingly, the agreement with the BCS results is quite 
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V. SUMMARY AND OUTLOOK 

In this paper, we aimed to describe the doping depen- 
dences of selected, experimentally measured quantities, 
we have decided to systematize the results coming from 
different versions of t-J model, that are discussed in the 
literature. From that analysis it follows that while the 
presence of the term ~ ci does not influence remarkably 
the results except for the largest doping x, the inclusion 
of the three-site terms {ci — C2 — 1) reduces substantially 
the range in which superconductivity is present. 

We have implemented the mean-field renormaliza- 
tion scheme^^ to t-J model, in which both the kinetic- 
exchange and the three-site terms can be taken into ac- 
count. Such RMFT approach, which is based on an ef- 
fective single-particle picture, is expected to be valid first 
of all in the overdoped regime, where an unconventional 
form of the Fermi liquid is obtained. The variational ap- 
proach based on the maximum entropy principle^^ has 
been used to determine mean-field parameters appear- 
ing in the model. The theoretical results yield a correct 
value for the upper critical concentration for the high- Tc 
d-wave superconductivity disappearance. 

Our method provides also lower single-particle energies 
compared to the previous MF results.— Consequently, 
a good estimate of the experimentally determined gap 
^k=(7r,o) is detected experimentally in the overdoped 
regime for typical values of the model parameters. How- 
ever, the values of the Fermi velocity in the nodal di- 
rection are still slightly too low. We also have exam- 
ined the temperature dependence of the superconduct- 
ing gap magnitude and have determined the critical tem- 
perature Tc{x) evolution as a function of hole concen- 
tration. As may be expected, mean-field results overes- 
timate Tc by a factor 3-5. The present results can be 
generalized by taking into account more complex lattice 
or band structure, the broken C^v symmetry (Pomer- 
anchuk instability^^), and antiferromagnetism. Further- 
more, study of the Fermi-surface-topology evolution with 
the doping is achievable. Finally, a more advanced 
scheme of calculating the entropy of the correlated state 
is desired to provide more satisfactory description of the 
nonzero-temperature situation. We should be able to see 
progress along these lines soon. 



FIG. 6: (Color online) Reduced temperature dependence 
of renormalized superconducting order parameter (A)c = 
Ac(r) for various doping levels. Squares: x = 0.125, solid 
circles: x — 0.25, diamonds: x = 0.3 are all for ci = 1, C2 = 0, 
t'/t = -0.25. Solid line: BCS result. 

good, despite the renormalization of A and its non- 
trivial d-wave symmetry. Uncorrelated (bare) gap ratio 
A(r)/A(0) exhibits very similar scaling with T/Tc- 
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